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Abstract. Let n particles move in standard Brownian motion in one dimension, 
with the process terminating if two particles collide. This is a specific case of 
Brownian motion constrained to stay inside a Weyl chamber; the Weyl group for 
this chamber is A„_i, the symmetric group. For any starting positions, we compute 
a determinant formula for the density function for the particles to be at specified 
positions at time t without having collided by time t. We show that the probability 
that there will be no collision up to time t is asymptotic to a constant multiple of 
£-n(n-i)/4 ag ^ g 0es £ infinity, and compute the constant as a polynomial of the 
starting positions. We have analogous results for the other classical Weyl groups; 
for example, the hyperoctahedral group B n gives a model of n independent particles 
with a wall at x — 0. 

We can define Brownian motion on a Lie algebra, viewing it as a vector space; the 
eigenvalues of a point in the Lie algebra correspond to a point in the Weyl chamber, 
giving a Brownian motion conditioned never to exit the chamber. If there are m 
roots in n dimensions, this shows that the radial part of the conditioned process 
is the same as the n + 2m-dimensional Bessel process. The conditioned process 
also gives physical models, generalizing Dyson's model for A n ~i corresponding to 
su n of n particles moving in a diffusion with a repelling force between two particles 
proportional to the inverse of the distance between them. 



1. Introduction 

Let n particles move in standard Brownian motion in one dimension, with the 
process terminating if two particles collide. Given the starting positions, we can 
use a reflection argument to calculate the density function for the particles to be at 
specified positions at time t without having collided by time t. Using this density 
function and the theory of Lie algebras, we can prove the following results; Theorems |3| 
and |] were first proved by Dyson [JTTJ . 



Theorem 1. For any starting positions, the probability that there will be no collision 
up to time t is asymptotic to a constant multiple of t~ n( ~ n ~ l ^l A as t goes to infinity; 
the constant is a known polynomial in the starting positions. 
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Theorem 2. Given that there is no collision up to time t, the distribution of the 
radius of the vector whose coordinates are the positions of the particles, divided by 
the square root oft, converges in measure to the distribution of the Bessel process with 
parameter n(n + l) / '2 at time 1, which is the radial part of an n{n + \) / 2- dimensional 
standard normal. 

Theorem 3. We can construct an n-dimensional Brownian motion which is condi- 
tioned for no two particles ever to collide. If we take Brownian motion on the space 
of Hermitian matrices, the induced process on the eigenvalues is the same process. 
If the starting point is appropriately chosen at the starting radius, the process given 
by the radial part of the conditioned Brownian motion is the same process as the 
n 2 -dimensional Bessel process. 



Theorem 4. This conditioned process is identical to the process obtained by n par- 
ticles moving in a one- dimensional diffusion with constant infinitesimal variance, 
with a repelling force between two particles proportional to the inverse of the dis- 
tance between them; that is, its infinitesimal generator has a 2 Ax,t) = 5ij, t) = 

This is a specific case of Brownian motion in a Weyl chamber; the vector whose 
coordinates are the locations of the n particles is constrained to stay inside the 
chamber. The Weyl groups is A n _i, the symmetric group. 

We have similar results for other Weyl groups; Theorems [I] and ||] have analogues for 
the classical Weyl groups, and Theorems |3] and f| have analogues for all Weyl groups. 
For example, the Weyl group B n models n particles in independent Brownian motion 
with an absorbing wall at x = 0. In this case, the asymptotic probability is ct~ n I 2 
(and the constant is again known), with the radial part corresponding to the n(n+l)- 
dimensional Bessel process; the conditioned process has a radial part which is the 
n(2n + l)-dimensional Bessel process. This generalizes the results of Pitman and 
Williams [p4 , 28j that one-dimensional Brownian motion conditioned to stay positive 
is the same process as the three-dimensional Bessel process. 

Our reflection argument is a generalization of the reflection principle, a standard 
argument in the analysis of both discrete random walks and Brownian motion. In 
the discrete case, it is used in the classical formula for the Catalan numbers, which 
enumerate the arrangements of n +l's and n — l's so that none of the partial sums 
are negative. Similarly, it can be used to study Brownian motion in one dimension 
with an absorbing barrier at x = and a known starting point . It has also been 



extended to multiple reflections in one dimension to study Brownian motion with two 
absorbing barriers [[□] . 

The reflection principle has been generalized to multiple dimensions. For example, 
the ballot problem, a classical problem in random walks, asks how many ways there 
are to walk from the origin to a point (Ai, . . . ,A n ), taking k unit-length steps in 



BROWNIAN MOTION IN A WEYL CHAMBER 



3 



the positive coordinate directions while staying in the region X\ > %2 > ■ ■ • > x n . 
The solution is known in terms of the hook-length formula for Young tableaux; a 
combinatorial proof, using a reflection argument, is given in [27|, [25]. 

The same reflection argument has also been applied to the case of n independent 
diffusions, or discrete processes which cannot pass each other without first colliding. 
Using this method, Karlin and McGregor jT7| give a determinant formula for the 
probability or measure for the n particles, starting at known positions, not to have 
collided up to time t and to be in given positions. Hobson and Werner [16J generalize 
this argument to n particles in an interval or circle, and use this to prove a result 
analogous to Theorem [| for n particles on the circle. 

Gessel and Zeilberger ||f 2|| , and independently Biane give a further generaliza- 
tion. For certain "reflectable" random walk-types, we can count the number of fc-step 
walks between two points of a lattice, staying within a chamber of a Weyl group, in 
terms of numbers of unconstrained walks. The steps must have certain allowable 



lengths and directions. In [13], all cases in which this method applies are enumer- 



ated, and determinant formulas are given for many important cases, including walks 
in the classical Weyl chambers. 

The argument of |12| can be generalized to Brownian motion in any Weyl chamber 
or chamber of a Coxeter group, with either absorbing or reflecting boundary condi- 
tions. We prove this generalization, and then use the result to compute determinant 
formulas for Brownian motion in the Weyl chambers of A n _i, B n = C n , and D n . 
The A n _i and B n cases are applicable to the independent motion of n particles in 
one dimension. The A n _i formula appears in [JPT| , using the model of n independent 
particles rather than motion in a Weyl chamber. The cases of the affine Weyl groups 



A n _i and B n are studied in |Ib| , also as models of n independent particles, in a circle 
or interval. 

These determinants factor into forms which can be easily analyzed; this allows us 
to find the asymptotic probability that there will be no collision up to time t for these 
three cases, for any starting point. In the case of A n _i, we get a simple formula for 
the actual measure as well. 

Weyl chambers arise naturally as the spaces of eigenvalues of elements of Lie al- 
gebras. Also, Brownian motion on a Lie algebra is symmetric under an orthonormal 
change of basis. There is thus a natural correspondence between standard Brown- 
ian motion on the Lie algebra and some diffusion on the Weyl chamber. Dyson flC 



showed that this diffusion has the same generator as the conditioned Brownian mo- 
tion for A n _i, this generalizes naturally to other Lie algebras. We can thus use the 
known properties of random matrices to study the distribution. 

DeBlassie || uses a different approach to give a more general formula for asymp- 
totics and density functions for a general class of cones, which include the cases 
discussed here. A cone is defined as the union of all rays from the origin which inter- 
sect the unit sphere in a connected open set C. The Laplace-Beltrami operator Lgn-i 
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on the unit sphere is the non-radial portion of the ordinary Laplacian. If Ai is the 
eigenvalue of Lgn-i with smallest absolute value on the space of all L 2 functions on 
C which vanish continuously on the boundary of C, then the asymptotic probability 
is a constant multiple of 

t - [-(n-2) + v /(n-2)2+4A?] /4 

This result gives asymptotics for a large class of cones. The coefficients depend on 
the eigenvalues of the Laplace-Beltrami operator rather than on explicit coefficients; 
these eigenvalues are known for Weyl chambers @. The asymptotics can thus be 
computed from these formulas; in theory, the explicit density functions can also be 
computed, but as infinite series, This formula also shows that the asymptotics for all 
such cones exist and are powers of t, with no other terms such as logarithms. 

The asymptotic probability is known to be a constant multiple of t~ m / 2 for a wedge 
of angle ir/m in two dimensions 0. The result holds in general, although the region 
is only a chamber of a Coxeter group (the dihedral group) if m is an integer. 

O'Connell and Unwin study the opposite asymptotic problem to ours, com- 
puting an asymptotic for the probability that n particles in independent Brownian 
motion will have a collision up to time t when t is small compared to the initial 
separation of the particles. 

Our results are organized as follows. Section 2 contains the basic definitions. In 
section 3, we prove the basic reflection result, and in section 4, we apply this result 
to get the determinant formulas. In section 5, we prove Theorems [l] and [|, and their 
analogues for B n and D n . In section 6, we construct the conditioned motion, both by 
Lie theory and by /i-transformation, proving Theorems || and |] and their analogues, 
and use the /^-transformation to find a physical model for all finite Coxeter groups. 

2. Definitions 

We will study a process with continuous sample paths in MJ 1 , either unconstrained 
or constrained by a chamber. In the constrained case, we may have either an absorb- 
ing boundary condition, causing the process to terminate when it hits a wall, or a 
reflecting boundary condition. All references to the analogous discrete problem are 
discussed in []13~| . 

We require that our chamber C be a chamber of a finite or affine Coxeter group. 
In the finite case, C is defined by a system of simple roots A C l n as 

(1) C = {f6l n |(a,f)>0 for all a G A}, 

and the orthogonal reflections r a : x t— > x— 2 ^^ a, generate a finite group W of linear 
transformations, the Coxeter group. In the affine case, the hyperplanes of reflection 
which define C do not all pass through the origin, and the group W is infinite, but if 
T is the subgroup of all translations, W/T must be finite. In the analogous discrete 
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problem, the steps of the the random walk must generate a lattice L which is stable 
under the action of W; in this case, C is a Weyl chamber and W a Weyl group. 

Let X(t) be a Markov process with continuous sample paths with values in R n ; that 
is, the distribution of X(t 2 ) given X(ti) is independent of X(t) at any £ outside the 
interval [£1^2] ■ We say that the constrained motion is reflectable if the increments 
of the unconstrained motion are symmetric under the Coxeter group; that is, the 
distribution of X(t 2 ) given X{t\) is the same as the distribution of w(X(£ 2 )) given 
the w-image w(X(ti)) as the starting point for any w G W. 

Standard Brownian motion is reflectable for any Coxeter group. For a finite Coxeter 
group, in which all planes of reflection pass through the origin, any diffusion with 
variance dependent only on time and the radius, and drift dependent on time and 
symmetric with respect to rotations and reflections about the origin, is reflectable; 
for example, there could be an absorbing or reflecting barrier at \x\ = R. 

As another example, consider the case in which each coordinate Xi(t) is an inde- 
pendent identical diffusion; this could model n independent particles instead of one 
particle in n dimensions. If our Coxeter group is the symmetric group A n —i — S n 
(giving the chamber x\ > x 2 > ■ ■ ■ > x n ), it permutes the particles, so the process is 
reflectable under this action; this case is discussed in [jT7|]. If the individual diffusions 
are symmetric about Xi = 0, then the product process will also be reflectable under 
the hyperoctahedral group B n , which includes all permutations with any number of 
sign changes. 

In the discrete case, reflectability requires the additional condition that the walk 



cannot go from inside the chamber to outside it without stopping on a wall ; this 
is our condition of continuous sample paths, which is satisfied by any diffusion. 

For fixed £, this process defines a probability measure Pt{A) = P{X(t) G A} which 
represents the chance that this process, if started at 0, will be in a set A at time 
£. We let c t (x) be the density function of this probability measure with respect to 
Lebesgue measure on M"; that is, 

c t (x)=P t (A). 

Now, we study the case in which the motion is constrained by a chamber, with ei- 
ther absorbing or reflecting boundary conditions. We must now fix the starting point 
77, since the process now depends on it in a non-trivial way. This generates new sto- 
chastic processes Y(j], t) for absorbing boundary conditions and Y'(jj, t) for reflecting 
boundary conditions. These give probability measures Q t (rj,A) and Q' t (r),A) which 
give the probability that the process, if started at rj, will be in the set A at time £, 
and these have density functions b t (rj, A) and b' t (r], A). Note that the total measure 
Qt{i]i will be less than 1 because the process terminates when it reaches a wall 
of the chamber. The total measure Q' t (r),M. n ) for reflecting boundary conditions will 
still be 1. 
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If our process is a Brownian motion, it has drift and variance of in each co- 
ordinate direction. If the Coxeter group contains any reflection in the Xj, x,,-plane 
other than a sign change of Xj or Xj, the reflect ability condition requires that <jj = o~j. 
Thus, for an irreducible Coxeter group, all of the (Tj must be equal; for a reducible 
Coxeter group acting on W 11 © IR™ 2 © • • • © M. Uk , we can multiply the coordinates in 
each R n J by a constant factor so that the <Tj are all equal. We can then re-scale time 
so that all Oi = \. 

Reflect ability also requires the Coxeter group to fix the vector // whose coordinates 
are the drifts /ij. In all non-trivial cases except for A n _\ on M. n , this requires that all 
Hi = 0, giving standard Brownian motion. For A n _i, it requires that the //j all be 
equal; we can then change coordinates to x\ = x, — /ijt to get an equivalent process 
in which all of the yJ i are zero. Thus, if the unconstrained process has stationary 



increments, and is thus a Brownian motion [|14j], we may assume that it is standard 
Brownian motion; however, we can state the theorems just as easily in terms of the 
more general reflectable process. 



3. The formulas 

Theorem 5. If Ct is the density function for a reflectable continuous stochastic pro- 
cess, then for absorbing boundary conditions, we have 

(2) b t (rj, A) = ^ sgn(»QO(A) - 77), 

and for reflecting boundary conditions, 



(3) K( V ,X) 



#{w e W : w(X) = A} 



If W is an affine group rather than a finite group, these may be infinite sums. The 
integrals over the images of any region must converge absolutely, since the measure 
of unconstrained motion over the whole space, the set of all IV-images of all points 
in the chamber, is 1. 



Proof. The discrete result analogous to (0) is proved in [12] and B; the proof which 



follows is essentially identical to that in [12] except that the discrete terms "walk" 



and "step" are replaced by their continuous analogues "path" and "time." 

Every path from rj to any w(X) which does touch at least one wall of the chamber 
has some first time to at which it touches a wall; let the wall be the hyperplane 
perpendicular to a«, choosing the largest i if there are several choices |2t|. Reflect 



the path after time to across that hyperplane; the resulting path is a path from r\ 
to r a .w(X) which also first touches wall i at time to- This clearly gives a measure- 
preserving bijection of paths, and since r a . has sign —1, all such paths cancel out 
in (□). The only paths which do not cancel in these pairs are the paths which stay 
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within the Weyl chamber, and since w(X) is inside the Weyl chamber only if w is the 
identity, this is the desired measure. 

For we note that the map on all paths starting at rj which takes every point 
to its unique image in the chamber C is measure-preserving, since we have reflecting 
boundary conditions and increments which are stable under the group W. This map 
takes all paths which end at any w(X) to paths which end at A itself. If A is on a wall 
of the chamber, paths to A may be counted multiple times, so we must divide by the 
size of its stabilizer in W. 

In practice, we can ignore this constant factor; it is 1 except for A on a wall, and 
this is a set of measure zero unless t = 0. Eliminating the denominator thus changes 
the density function only on a set of measure zero, and thus does not change the 
measure of any measurable set. 

4. Determinant formulas for the density functions 

We can now apply this theorem to standard Brownian motion, in the Weyl cham- 
bers of A n _i, B n = C n , and D n , with either absorbing or reflecting boundary 
conditions. The measure for unconstrained standard Brownian motion is = 
Jir=i Nt(xi), where N t is the normal distribution function with mean and variance 
t. Since this factors into separate terms for the individual coordinates, we can use 
the same techniques to compute determinant formulas as in the discrete case |I3|. 



The most interesting case is A n ^\ = S n , the symmetric group. The Weyl chamber 
is x\ > X2 > ■ • ■ > x n . This Brownian motion thus models n independent parti- 
cles in one dimension. With absorbing boundary conditions, collisions are forbidden 
(the process terminates if one occurs); with reflecting boundary conditions, particles 
collide elastically with one another. 

For absorbing boundary conditions, we write the sum (0) as 

(4) b t (7], A) = J2 sgn(<7)ct(<7(A) - 77) 
and use the value of c t to write this as 

n 

(5) b t ( V , \)=Y< II ( N t(^d) - *)) ■ 

<T(zS n 1=1 

This sum can be written as a determinant, which gives 

(6) 6 t (7 ? ,A) = det|JV t (A i -77 J -)|. 

nxn 

This determinant gives the measure for n particles which start at positions rji and 
are in independent Brownian motion to be at positions Aj at time t without having 
collided. 

For a more general product of n independent diffusions, with individual density 
p t (x — > y) for the diffusion which started at point x to be at point y at time t, we 
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can apply the same argument. This gives the following generalization of 



which 



first appears in [T7 



(7) b t ( V ,X) = det\p t (r ]j ^X i )\. 

nxn 

For reflecting boundary conditions, the calculations are the same except that there 
is no term for sgn(cr), and thus we get permanents rather than determinants. This 
result can also be seen by observing that an elastic collision between two identical 
particles is equivalent to the two particles passing through each other with no collision, 
and thus the reflected particles will be at the positions indicated by A if the unreflected 
particles are at any permutation of the coordinates of A. 

For B n , the hyperoctahedral group, which includes permutations with any number 
of sign changes, the Weyl chamber is X\ > £2 > • • • > x n > 0. This also models n 
independent particles in one dimension, with an additional wall at x — 0. 

We write w G W as a product of an e which negates some coordinates and a a in 
the symmetric group. We get 

n n 

(8) b t ( V ,\)(x)=Yi e sgn^n^n^^w-^))- 

a&Sn e<=±l i=l i=l 

Using the multilinearity of the products in the determinant, we can again write this 
sum as a determinant, with separate terms for e« = 1 and q = — 1 in each entry. We 
use N t (x) = N t (—x) to keep the signs of A positive and get a more elegant formula. 
This gives 

(9) bt(ij, A) = det \N t (Xi - rjj) - N t {Xi + Vj )\ . 

nxn 

This determinant gives the measure for n particles which start at 77, to be at Aj at 
time t, neither having collided nor having touched x — 0. 

Again, the same argument applies if n particles are in general independent diffu- 
sions, provided that the diffusions are symmetric about x = 0. The more general 
formula is 



(10) b t {r], A) = det \p t (r)j -> Aj) - p t (rjj -> -Aj)| . 

nxn 

For reflecting boundary conditions, we lose the sign of the a, which makes the 
determinant into a permanent, and the sign of the £j, which turns the minus sign 
between the two N t in (||) or p t in (|T0"|) into a plus sign. Again, the resulting formula 
is the same as would be obtained by treating elastic collisions as though the particles 
passed through each other, and allowing particles to pass through the wall at x = 
instead of bouncing. (In the transformed model, particles at positions x and — x no 
longer collide, but since they collided elastically in the original model and could be 
considered to pass through each other instead, the effect is the same.) 
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For D n , the even hyperoctahedral group, which includes permutations with an even 
number of sign changes, the Weyl chamber is x\ > x 2 > ■ ■ ■ > x n , > — x n . This 
does not give a natural model for n particles in one dimension. 

Again, we can write w = ea. We take our sum over all possible e, and then add 
an additional factor of (1 + Yl™ =1 q)/2 to annihilate those e which are not allowed in 



(11) 



1 _|_ FT" n 

b t ( V ,X)(x) =J2J2 sgn((7) + 1 6i J] (uNtfaKa - Vi}) 



crGSn £i=±l 



1=1 



We now take the | and the dllLi e «)/2 terms separately. The dllLi e «)/2 term is half 
the sum we had in (||); the 1 term gives half of (||), but with a plus sign between the 
terms. Thus we get 



b t (rj, X)(x) 



1 r 

2 



(12) 



det \N t (\i - rjj) - N t (Xi + rjj)\ 

nxn 

+ det iNtiXi-rj^+NtiXi+rjj)] 



If we let A' be obtained from A by changing the sign of X n , this will change the sign of 
the first term but preserve the second term. Thus the first term alone, with no factor 
of 1/2, is b t (r], X)(x)—b t (r], A') (a;), and the second term alone is b t (r], X)(x)+b t (r], X')(x). 
If A n = or r] n = 0, the first term is zero, so b t (i],X)(x) is the second term alone, 
with the factor of 1/2. 

For reflecting boundary conditions, we ignore the sign of a. Thus the determinants 
become permanents, but we keep the minus signs because they came from the factor 
(1 + nr=i e «)/2 ; which was not from sgn(w). 



5. ASYMPTOTICS 

We can use these formulas to find asymptotics for the probability that the motion 
will not hit a wall of the chamber by time t, and for its distribution at time t given 
that it has not hit a wall. 



5.1. Calculating the individual values for A n -±. We can eliminate the determi- 
nant to get a more explicit formula for (|6|) at a single point A if all the coordinates of 
our starting point r\ are rational. Re-scaling by t — > c 2 t will make all the co- 

ordinates integers, and we can then translate all coordinates by — r] n so that we have 
i] n = 0. (Both of these transformations leave the Weyl chamber x\ > x 2 > • • • > x n 
unchanged.) 
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We now write out the normal distributions in (0) explicitly as exponentials, and 



expand (Aj — rjj) 2 as Af — 2XiT)j + rj 2 : 



(13) 



- : — 777 det 

(27Tt)"/ 2 nxn 



exp 



-A 2 + 2X i 7 lj - j]) 
2t 



Row i of this matrix contains a constant factor exp(— X 2 ), and column j a constant 
factor exp(— rj 2 ), so we can take these out, and put them in a constant term, which 
simplifies further because |A| 2 = ^ Af . This gives us 



(14) 



(27rt)"/ 2 



exp 



IAI 



2t 



det 



exp 



XiJJl 
t 



Since the rjj are all integers, we can write the determinant as the generalized 
Vandermonde determinant 



(15) 



det |[exp(V*)] 



If rjj — n — j, this is the standard Vandermonde determinant, equal to 
(16) IJjexptVO-exp^./O]- 



i>j 



And for any non-negative integers r]j, it is the product of this Vandermonde deter- 

,»?n(exp(Ai/t),... ,exp(A n /t)) 



minant and the Schur function 



s, 



m- 



-n+1,772— n+2,. 

The Schur function can also be defined combinatorially |I9| , with the coefficient of 
Y[ xT in Sfi the number of ways to fill in the partition diagram of fi, using the number 
i exactly rii times, such that the entries are non- decreasing in each row and strictly 
increasing in each column. 

In particular, we can let C v = S a (l, 1, . . . , 1) be the total number of such tableaux; 
this is important because it is the approximate value of the Schur function when the 
Aj are much less than t. (This holds because the Schur function is a homogeneous 
polynomial of degree t]i + ■ ■ ■ + i] n — (n(n — l))/2, with positive coefficients.) This 
will show the dependence of the asymptotics on the starting point. 
This constant is known 



19 ; it is 



(it) a^n^-^/nu-z). 

i<j ' i<j 

This allows us to compute the constant term in the asymptotics. 

For any 7] with the same sum of the coordinates, all multiples of 1/c, our rescaling 
gives a Schur function whose index is the transformed vector a]. Rescaling to restore 
the old time values gives a Schur function which is a homogeneous polynomial in the 
exp(Aj/t) 1//c whose degree in the exp(Xi/t) is 771 + • • • + r] n — [n{n — l))/2c. Thus, 
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for any point 77, even one with fractional coordinates, the degree will be bounded by 
the sum of its coordinates. Thus the determinants, and therefore the b t (r),X), for 
different starting points rj and rj' both of radius less than a known 5, will be in the 
approximate ratio of C v to Cy, with an error of 0(8/ 1). If 5 is fixed, then when we 
calculate the asymptotics as t — > 00, the Schur function will converge to C v at a rate 
ofO(l/t). 



5.2. Asymptotic probability of no collisions. The integral of the value in (|b 
over the whole Weyl chamber is the probability that n particles starting at the po- 
sitions r]j will have no collisions up to time t. We can use this formula to show 
that the asymptotic probability as t — > 00 is a constant multiple of £- n (™- 1 )/ 4 ) with 
the constant depending on rj. We will also show that, given that the n-dimensional 
Brownian motion has not hit a wall, its radial distribution, rescaled by multiplying 
the radius by 1/ converges to the distribution of the Bessel process with parameter 
n(n + l)/2 at time 1. 

Fix t very large compared with \i]\ 2 . If | A | 2 is much larger than t, then the ex- 
ponential of — \X\ 2 /2t in (14) will decay exponentially fast. In particular, \/t can 
be assumed to be O^ 6-1 / 2 ), so we only need to look at the leading nonzero terms 
in fll5|); the error will be a factor of this order when compared with the value of 
the determinant. In (pf), the first nonzero term of the factor [exp(Aj/t) — exp(Xj/t)] 
is (Aj — Xj)/t; we can then take out the factor of f-^- 1 )/ 2 an d leave only a term 
involving the A. Likewise, in the Schur function from fll5|), we need only keep the 
constant term C v . 

Thus our probability is asymptotic to the integral over the Weyl chamber of 



(2-nt) n /' 1 



(- 


X 


2 _ 


V 


to 


\ 2t 



i>j 



And this integral can be computed by using Selberg's integral 20, 22]; we have 



dxi = 2 3 ™/ 2 JJr((fc/2) + l). 

k=l 

This corresponds to our desired integral when we set x = X/y/t, and divide by n\ 
because we are taking our integral over only one of the n! different Weyl chambers. 
Dropping the exp(— \r)\ 2 /2t) (which goes to 1 as t — > 00 and thus doesn't affect the 
leading term), and writing out C v explicitly again gives our full asymptotic: 

(20) IL<jfa-^) (27r) n/2 2 3n/2 JJ[r((ifc/2) + l)}^^^. 



(19) 



exp(— |x| 2 /2) 



Uf 

i>3 
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We can also note that the exponential in ( jl8[) is spherically symmetric, while 
Ili<j(Aj — Aj) is homogeneous of degree n(n — 1)/2. Thus the density, integrated over 
the sphere of radius r at a fixed t, is a constant multiple of 

exp(-r 2 /2t)r«"("- 1 )/ 2 )- 1 , 

and thus a constant multiple of the radial distribution of (n(n + l)/2)-dimensional 
Brownian motion at time t. Thus, given that no two particles have collided, the 
distribution would be exactly the same if it were given by ([18]). This proves Theo- 
rem § if we restrict to |A| < ct 1 ^ 2+e , the \/t terms are all 0(t e ~ l l 2 ). Thus, in this 
region, the ratio of the radial distribution for the constrained motion to the radial 
distribution for unconstrained motion in n(n + l)/2 dimensions converges uniformly 
to 1 at a rate of 0{t e ~ 1 / 2 ). 

Equivalently, we could fix time and A, and for a scalar 5, take rj = Srj as our 
starting point. By the scaling properties of Brownian motion, this is equivalent to 
keeping rj fixed, taking A = X/5, and t = t/5 2 . Thus, as 5 — > 0, the ratio of the 
radial distributions converges uniformly to 1 within the region |A| < ct 1 ^ 2+t at a rate 
of 0(5 l ~ € ), and the probability that either distribution is outside that region goes to 
exponentially fast. 

5.3. Asymptotics for B n : no collisions and a wall. We can use the same tech- 
nique to get asymptotics for B n as for A n _i, here, we get a constant multiple of t~ n I 2 
as the probability of no collision, and the radial distribution given no collision is the 
n 2 + n-dimensional Bessel process. 

Here, we require that the coordinates all be odd integers. Again, we write out the 
determinant (0) explicitly: 



(21) Un. A) = r^rdet 



nxn 



-A? + 2A^ -rj]\ ( -A? - 2A^ - rj 



exp — exp 



2t J V 2t 



As before, we remove the constant factors to get 



(22) b t ( V , A) = j—^j- 2 exp ( L_ ) d et 



(- 


A 


2 _ 


V 


\ 2t 



nxn 



exp l — exp ' 



The determinant here is not an actual Vandermonde determinant. However, in 
the specific case rjj = 2n + 1 — 2j, we can make it a Vandermonde determinant by 
elementary operations. Adding (— l) fc ( 2n+ ^~ 2j ) times column j + k to column j does 
not change the determinant, but it changes the entries in column j to 



2rt+l-2j 



Zz~ J /o _|_ i _ o '\ 

E ("!)* ( n k 3 ) ex P(( 2 ^ + 1 - 2J - 2*) V*) 

i — n \ / 

= [exp(A i A)-exp(-A l /t)] 2ri+1 ^'. 
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This is a generalized Vandermonde determinant; we can make it an actual Van- 
dermonde determinant by dividing row i by exp(Aj/i) — exp(— Xi/t). The resulting 
determinant is 

(23) det | [exp( A*/*) - exp^/i)] 2 ^ 1 ^ | , 

nxn 1 1 

and its value is 

(24) J] [(exp(Vt) - exp(-A,A)) 2 - (exp(A J /t) - exp(-A,/t)) 2 ] . 

i>j 

Putting this together with the constants we have taken out, we get 
b t(v,*) = ,*ws exp 



(27Tt)"/ 2 ^ V 2t 

n 

JJ[exp(Ai/t) - exp(-V*)] 
i=i , 

II [( ex P( A *A) - exp(-A,/t)) 2 - (exp(A i /t) - exp(-A J /t)) 2 ] 



x 

i>j 

For a more general starting point, with all coordinates odd integers, we note that 
(exp(\iT)j / 1) — exp(— Xi7]j/t)) is a polynomial in exp(Aj/t) — exp(— Xi/t) with no con- 
stant term. We can break these polynomials into their individual terms, giving a 
large number of determinants, each one a generalized Vandermonde determinant in 
the exp(Aj/t) — exp(— Xi/t). Each individual determinant is thus the product of fl24|) 
and a Schur function; it also contains the product of the exp(A«/t) — exp(— Xi/t) as a 
factor, since these are constant factors in row i. Thus each determinant is a product 
of these factors with some symmetric function in the exp(Aj/i) — exp(— Xi/t). As with 
A n _i, the error we get in approximating the Schur functions by their constant term 
C v is a factor of 0(l/t). 

And as with A n _i, we can get asymptotics by integrating this over the Weyl cham- 
ber. For large t, the radial exponential will be exponentially small if |A| > t 1 ^ 2+t , so 
we can assume that all of the Xi/t are very small. Thus exp(Aj/t) — exp(— Xi/t) can 
be approximated by its leading nonzero term, 2Xi/t. Thus, as with ($M, we get 



1 / I \ 1 2 I 1 2 \ n 

(25) (2^7i (^ir^ 1 ) w n w n<^ - *s 
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This integral can also be computed by using Selberg's integral POL p2| ; we have 



exp(-|x| 2 /2)J] 



r 2 _ T 2 ' 



o(rt 2 +3n)/2 n 

n[ r ( 1 + V2)r((i + *)/2)]. 



7T 



ra/2 



i=l 



This corresponds to our desired integral when we set x — X/yt, and divide by 2 n n\ 
because that is the number of Weyl chambers. Dropping the exp(— \r]\ 2 /2t) because 
it goes to 1 as 1/t, we get our asymptotic probability that Brownian motion started 
at the specific point rj will not hit a wall up to time t: 

n(3n 2 /2) n 

(26) J] [r(l + k/2)V{{\ + k)/2)\ t 



-n 2 /2 



TC' L n\ 



k=l 



For a more general starting point, each individual determinant in the sum gives an 
integral of a homogeneous polynomial in the exp(A,/t) — exp(— Xi/t) which is of degree 
at least n 2 , since it contains the previous determinant as a factor. The same technique 
as above gives an asymptotic which is thus at most t~ n I 2 . Thus the asymptotic is at 
most a constant multiple of t~ n / 2 , and is less if and only if the coefficient of t~ n I 2 
is zero. But the coefficient cannot be zero. Let m = max((r]i — i]i + i/2),i] n ), and 
ri'i = (2n + 1 — 2i)m. Then we know that Brownian motion starting at rf decays 
asymptotically as t~ n I 2 , but Brownian motion starting at rj will always hit a wall if 
it is translated to start at rj', because rj 1 is at least as far as rj from every wall. Thus 
the asymptotic probability will be some constant C v times the formula (|26|); we will 



compute C v in Section |5\5| . 

A result analogous to Theorem |] also holds for B n . The exponential in numbb- 
nasymp is spherically symmetric, while the product is homogeneous of degree n 2 . 
The same argument as for A n _\ shows that the ratio between the radial part of 
this distribution and the radial distribution for unconstrained Brownian motion in 
n 2 + n dimensions converges uniformly to 1 at a rate of 0(t €_1//2 ) inside the region 
|A| < ct 1//2+€ , and the probability that either distribution is outside that region goes 
to exponentially fast. Again, we can rescale by multiplying the radius by l/y/i to 
get convergence to a fixed distribution. 

5.4. Asymptotics for D n . The process for D n is almost the same as for B n , so we 
won't work it out in full detail; we get a similar result with the same error terms and 
convergence properties. 

Here, it is most natural to let rji = n—i. For this value of r], we have only the second 
determinant in (|1^), with a plus sign between the terms; the last row of the other 
determinant is zero. For general r], the first determinant is the B n determinant, which 



BROWNIAN MOTION IN A WEYL CHAMBER 



15 



we know is 0(t~ n2 / 2 ), and we will show that the second determinant is asymptotically 
larger. 

The determinant that we get is 



(27) 



det 



cxp 



+ exp 



Again, this isn't a Vandermonde determinant, but exp(Aj/t) +exp(— \i/t) n is equal 
to exp(\ir]j/t) + exp(—XiT]j/t) plus a sum of lower order terms, so elementary opera- 
tions which do not change the determinant give us the Vandermonde determinant 

(28) det I [exp(Ai/t) + exp(-A i /*)]"~" i l , 

nxn 1 1 

and its value is 

(29) ]J [exp(Vt) + exp(-A i /i) - exp(A i /t) - exp(-A,/t)] . 

i>j 



The leading nonzero term is 



i>j 



and the full integral is 



(30) 



(27Tt)"/ 2 



exp 



■|AP 



2t 



(2A)" 2 -"JI(A- -A|). 



i>j 



Again, we get a result which can be obtained from Selberg's integral p(| |22 
have 



wc 



(31) 



exp(-|a;| 2 /2) 



n 

i>j 



2(n 2 +2n)/2 



7T 



n/2 



n[r(i + i/2)r(i/2)] 



i=i 



(32) 



Here, there are 2 n 1 n\ Weyl chambers, which gives the asymptotic 

2(3n 2 -3n+2)/2 



X\[T{l+k/2)T{{l + k)/2)\t 



(-n 2 +n)/2 



fe=l 



For a general starting point whose coordinates are all integers, we use the same 
technique as for B n . The terms in the determinant are all polynomials in exp(Aj/t) + 
exp(— Xi/t), so we can again split the sum into individual determinants, each of which 



is the product of (|29|) and a Schur function of the exp(Aj/t) + exp(— Xi/t). Since Xi/t 
is small, these are all close to 2, so the full asymptotic is a sum of terms of order 
£(-n +n)/2_ ^ g w ifo -Q n ^ we can translate to i]' which is further from any wall than 7] 
to show that the coefficient of t <y ~ n +n )/ 2 cannot be zero. 
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5.5. The constant factor for a general starting point. Since we know the con- 
stant factor in the asymptotic probability of no collision up to time t for one specific 
starting point r), and the asymptotic distribution for an arbitrary starting point, we 
can use the time-reversibility of Brownian motion to compute the asymptotic prob- 
ability of no collisions for an arbitrary starting point. The argument is the same for 
all of the Weyl groups. 

In each case, the density b t (r], A) for Brownian motion starting at rj to be at A at 
time t, not having collided with a wall up to time t, is asymptotic for large t to a 
product of the form 

-iai 2 - \vr 



(33) C n C'J(t)exp 



2t 



in which C' x and f(t) are known. Since Brownian motion is symmetric in time, the 
density b t (\,rj) must be equal to b t (r], A). As long as t is large enough compared to 
| r) | and |A| for the formula (|33| ) to be valid (which it will be for large t because of the 
exponentials), we can reverse the roles of r\ and A. Thus we have C V C' X = C\C' n) and 
since the formula for C is known, we see that C v must be a constant multiple of C ; 
it is thus a constant multiple of YliVi ~ Vj) f° r An of Ylivf ~ v]) YlVi f° r B n) and 
Ilfo? - rg) for D n . 

We know the value of the constant from the formulas (|20|), (f26|), and (0). For 
A n , we already have the value for a general starting point because we used the Schur 
functions to compute C v in (p0|); we could have instead used this technique. For B n , 

we get 

2 2 n 

n ^ n 

11 (On - 1V2 _ O.i - 1^2 11 



i<j 
2 (3n 2 /2) 



(2j - l) 2 - (2z - l) 2 i A 2n + 1 - 2i 

n 2 /2 



^-^r\\w + k/2)m + k)/2)}t 

k=l 



as the asymptotic probability of no collision; for D n , we get 



n 



. (2j - l) 2 - (2i - l) 2 

o(3n 2 -3n+2)/2 n 

x ^ny^ + fc/2)r((i + ao/2)] t^ 2 +")/ 2 . 

fe=l 

6. Random matrices and conditioned Brownian motion 

6.1. Brownian motion on a Lie algebra. Instead of viewing our Weyl chamber 
as the chamber of a Weyl group, we can view it as the space of eigenvalues of the Lie 
algebra corresponding to that Weyl group, and then use the theory of Lie algebras 
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to study it. The following construction was first developed by Dyson [|T0| for su n ; he 
computed the properties of the Brownian motion from the specific data rather than 
using the general Lie theory. The results we use from Lie theory are given in p), 



Given a Lie group G and its Lie algebra g, we can define a normal distribution or 
Brownian motion on g by viewing it as a vector space and taking an orthonormal 
basis with respect to a G-invariant inner product. In particular, on so n (M), the Lie 
algebra of skew-symmetric matrices, we can take l/y/2 times standard Brownian 
motion on each matrix entry My with i < j, and take Mji = —My-; this makes the 
Hilbert-Schmidt norm of M equivalent to the radius of a Brownian motion in n 2 — n 
dimensions. Likewise, for su n (C), the Lie algebra of skew-Hermitian matrices, we can 
take standard Brownian motion on the imaginary part of each diagonal entry Ma, 
and l/y/2 times standard Brownian motion on each matrix entry My with i < j, 
with Mji = —Mij] this again gives us the Hilbert-Schmidt norm, this time in n 2 
dimensions. For any Lie algebra, the Brownian motion at time t will have a Gaussian 
distribution on the space of matrices in that Lie algebra. This allows us to use all of 



the known results about random matrices p2| to study the motion, and its eigenvalues 
in particular. 

Since this Brownian motion is invariant under conjugation by G, it induces a 
diffusion on the space of eigenvalues of the matrices. The Weyl chamber is in a 
correspondence with the space of eigenvalues, obtained by dividing the independent 
eigenvalues by i and arranging them in decreasing order. To study this diffusion, 
we consider the measure on the Weyl chamber which is induced by the standard 
measure on the Lie algebra. We can study this measure by proving a version of the 
Weyl Integration Formula for Lie algebras, which can be proved directly from the 
analogous formula for Lie groups . 

Lemma 1. Let g be a finite- dimensional Lie algebra with m roots, and f{X) = 
f(X\, . . . , An) a function on g with compact support which depends only on the eigen- 
values of X. Let 5' be (27r) m times the product of all the roots of g. Then the integral 
of f over g is equal to the integral ofd'S'f over the Weyl chamber. 

Proof. Let G be the Lie group corresponding to g. We can apply the exponential 
map h : X — > exp(iritX) to translate a neighborhood of zero in g which contains the 
support of / to an arbitrarily small neighborhood of the origin in G by choosing t 
sufficiently small. The derivative of this map can be made arbitrarily close to (nit) 11 
(which becomes (vrt) n when we consider volume elements) by taking t sufficiently 
small. We now apply the Weyl Integration Formula to h o /, which is a class function 
on the Lie group G; it says that the integral of h o / over G is equal to the integral 
of 55 over a maximal torus, where 

(34) 5 = Y[ (e™ (a?) - e -™ w ) . 
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For x sufficiently close to zero, we can approximate 5 by its leading nonzero term, 
which is 



(35) 5^5' = (2n) m JJ 



a{x). 

aGA 

This allows us to proceed with the integration. 

We can apply this formula to the characteristic functions of small intervals for 
the eigenvalues; this shows that the measure at a point (xi, . . . ,x n ) in the Weyl 
chamber is proportional to 55. In particular, we note that 5 is positive real on the 
interior of the Weyl chamber, and zero on the walls; thus we have a diffusion which 
is conditioned never to leave the Weyl chamber. 

Note that this proves Theorem |^ for this process, and an analogous theorem for 
any Weyl group. If we view g as a vector space, the Hilbert-Schmidt norm of a 
matrix is the radius of a vector. If there are m positive roots, then g has dimension 
n + 2m [|15||. The Hilbert-Schmidt norm is thus given by the Bessel process in n + 2m 



dimensions. The Hilbert-Schmidt norm is unchanged under diagonalization, so the 
diffusion in the Weyl chamber has its radius given by the same process. 

To properly state this result for Brownian motion conditioned never to leave the 
Weyl chamber, we cannot allow the process to start at the origin, which is not a point 
in the chamber. However, the square of the product of the roots is a homogeneous 
function of degree 2m; it thus gives an identical distribution on any fixed radius. 
Thus, if we start the motion on the Lie algebra at 0, then at any later time, given the 
fixed radius, the distribution on the sphere will be given by the square of the product 
of the roots, and the process can be continued from that time on as a Bessel process. 
This allows us to state the general theorem, formalizing Theorem |3]. 

Theorem 6. For any Weyl group W acting on W 1 with m roots, consider the process 
which starts at a fixed radius r Q , with the starting point chosen on the sphere of radius 
ro by a distribution with density proportional to the square of the product of the roots. 
Then the radius of the position of transformed motion at time t gives a Bessel process 
with parameter n + 2m. 

The simplest case of this theorem is for B\, which is one-dimensional motion with 
an absorbing boundary at x = 0. In this case, the only root is x, so m — 1, and we 
have the result of and that Brownian motion conditioned never to hit is the 



same as the three-dimensional Bessel process. In our more general cases, the number 
of dimensions for the Bessel process is n + 2m = n 2 for A n _i, n(2n + 1) for B n , and 
n(2n — 1) for D n . 

6.2. Construction by /^-transformation. We could have alternatively constructed 
this diffusion by /^-transformation ||. We will first construct the conditioned Brown- 
ian motion in this way, and then show that the two processes are actually identical; 
this allows us to develop the physical model of Theorem 01 
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6.2.1. General properties. Given standard Brownian motion in any number of dimen- 
sions, we can use the process of /i-transformation to construct a Brownian motion 
satisfying certain conditions. For any non-negative harmonic function h, the measure 
for the transformed Brownian motion to go to A at time t after starting at rj is h(X) 
times the measure for untransformed Brownian motion with the same starting point, 
renormalized by an appropriate constant so that the total measure on all paths is 1. 

In particular, suppose that we have a connected region D and a harmonic function 
h which is zero on the boundary of the region and positive on the interior. Then 
the /i-transformed Brownian motion will be conditioned to stay in the interior of the 
region. 

The density function for this transformed motion will be h times the density func- 
tion for the untransformed motion on the same region, normalized appropriately. 
That is, its value at a point A and time t will be h(X) times the value b t (r), A) which 
gives the measure for unconstrained motion to go from rj to A in time t while staying 
within the region, since the transformation puts a new measure on the same set of 
paths. 

The same technique can also be applied to a discrete random walk with a set S of 
steps, provided that h is harmonic in the discrete lattice; that is, we need 



The transformed discrete walk now has probability h(x+s)/(h(x)\S\) instead of l/\S\ 
of going from x to x+s in a given step. Thus, given a starting point rj, the probability 
of going to A in a given number of steps of the transformed random walk is h(\)/h(r)) 
times the probability for the untransformed walk to go to A while staying within the 
region in which h is positive. 

In order to make it impossible for the walk to leave the region D, we need h to be 
zero on all points which can be reached from the interior of D in a single step. Thus, 
only if the discrete walk is reflectable (as defined above) is it sufficient for h to be 
zero on the continuous boundary of D. 

6.2.2. Finding the function. The properties of conditioning, as well as our asymp- 
totics, suggest that 5 itself should be our h. It can be checked algebraically that h is 
harmonic for each group; however, it can also be proved naturally. 

Theorem 7. For any finite Coxeter group W , the product of all the positive roots is 
a harmonic function, both for the continuous Laplacian and for the discrete Laplacian 



(36) 





for any set S which is symmetric under the group W. 
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Proof. This result in the continuous case is due to 0; the discrete argument is a 
simple generalization which is mentioned in |J. 

By the properties of root systems, a reflection in any simple root changes the sign 
of only that root, while permuting the other positive roots. Thus the product h of all 
m roots is antisymmetric in all the simple roots, and by applying repeated reflections, 
we see that it is antisymmetric in every root. It is also of degree m. 

The continuous Laplacian is spherically symmetric, and thus symmetric under W . 
The discrete Laplacian is symmetric under W because the set S is. Applying the 
continuous Laplacian to a polynomial decreases the degree by 2, while applying the 
discrete Laplacian decreases the degree by at least 1 since h(x + s) — h(x) is of lower 
degree than h. Thus the application of the Laplacian to h gives a polynomial which 
is of degree less than m which is still antisymmetric in W. 

Now, any polynomial which is antisymmetric in W must be zero on every one of 
the hyperplanes of reflection. If it is not identically zero, it must have all m roots as 
factors, so it must be of degree at least m. Thus the Laplacian must annihilate our 
polynomial h, so h is harmonic for either the discrete or continuous walk. 

It follows that the transformation by this h gives the same process as the process 
generated by Lie theory. We have already computed the asymptotic density that 
Brownian motion started at a fixed rj will remain in the Weyl chamber for time t and 
be at A at that time. Transforming by h has the effect of multiplying the measure of 
all paths from a fixed 77 to an arbitrary A which stay within the chamber by a factor 
of h(X) (and a normalizing constant). Thus the density function for the transformed 
Brownian motion which starts at rj to be at A at time t is proportional to h(\)bt(r), A), 
and this converges to a constant multiple of h(X) 2 as t becomes large. This is the 
same factor S'S' which we obtained from the Weyl Integration Formula. 

6.3. The infinitesimal generator and physical models. We have constructed 
the conditioned process as a transformation of Brownian motion. We can also con- 
struct it as a diffusion with its known infinitesimal generator. In this form, both the 
multidimensional models and the models of n independent particles lead to natural 
physical models. 



We will use the notation of |18| for infinitesimal generators of diffusions. The drift 
vector n is defined by 

t) = hm ^{E(Xi(t + At)) - Xi(t) \X(t) = x}, 

and the infinitesimal variance matrix is defined by 

njfrt) = hm ^{[E(Xi(t + At)) - XmiEiX^t + At)) - X 3 {t)]\X{t) = *}, 

We will omit the variable t in the infinitesimal generators of our diffusions, because 
they are independent of time. 
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It can easily be checked that Brownian motion transformed by a harmonic (and 
thus necessarily C 2 ) function h has infinitesimal drift // = Vh/h and infinitesimal 



variance = 5y. 

Note that \7h/h = V(log/i). This makes the computation easy, because our h is 
the product riaeA a °f a ^ the ro °t s ; viewed as linear functions of the X{. If we write 
a(x) as the dot-product (a,x). then we have 

(37) = ^ V (a,a) = 5: W 

aeA agA v ' ' 

Since (a,af)/|a| is the distance from x to the hyperplane orthogonal to a, while 
a/\a\ is the unit vector in the direction of a, this term in the drift is the inverse of 
the distance between x and the hyperplane, directed away from the hyperplane. We 
thus have the following physical model. 

Theorem 8. For any finite Coxeter group, Brownian motion conditioned to stay 
within a chamber is equivalent to the motion of a particle in a diffusion with constant 
infinitesimal velocity, and a repulsive force from every hyperplane of reflection (not 
merely the walls of the chamber) inversely proportional to the distance from that 
hyperplane. 

We can also look at the n coordinates as individual motions in one dimension. If 
the root a contains cxj, we get a term c/(a,x) in the sum. 

In particular, for A n , the drift /ij is Yu&j ~ x j)- Thus each particle is subject 
to a repelling force from every other particle (not merely its neighbors), inversely 
proportional to the distance between them. This proves Theorem [|; this result is 
originally due to Dyson [T^ 



The model for B n is not as natural as a model of particles. The drift /ij is 

1 



+ 



ry . V y» . rr> . ry . I rt~> 



#3 

That is, each particle is repelled by every other particle, and by the wall at (the 
1/xi term), but also by the mirror image of every other particle (the l/(xi+Xj) term), 
as if the wall at x = was also a mirror reflecting all forces. For D n , we have only 
the terms of l/(sc» — Xj) and l/(xj + Xj); this means that the mirror reflects forces 
but is itself permeable to particles. 

These models are more natural as models of the eigenvalues. For D n , a matrix in 
the Lie algebra S02 n has eigenvalues ±i\j] this model thus says that the eigenvalues 
in different pairs ±iAj and ±i\k for j ^ I repel one another other, although the 
pair ±iAj do not repel each other. For B n , a matrix in the Lie algebra S02n+i has 
eigenvalues ±i\j and 0; this model thus says that the unconnected eigenvalues i\j 
and i\k for j ^ k repel each other, and each eigenvalue is also repelled by the fixed 
eigenvalue at 0. 
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7. Open problems 

The discussion of Brownian motion on a Lie algebra is valid for the exceptional 
Lie groups as well, but the techniques for computing the specific asymptotics do not 
appear to work. For a chamber of a general Coxeter group, the Lie algebra technique 
is not meaningful. In either case, is it possible to get the same type of asymptotics, 
with constant terms in particular? 

We have shown that the same harmonic functions which we used to transform 
Brownian motion can be used to transform discrete random walks. Is it possible 
to use these results to compute asymptotics for the discrete walks, including the 
constants on leading terms? 



Brownian motion can be defined on a general manifold • This allows us to apply 
the argument of Theorem [| whenever we have a suitable chamber. As before, the 
Brownian motion must be symmetric under the reflections in any wall of the cham- 
ber, and the reflections in the walls must generate a discrete group which partitions 
the manifold into chambers. For example, since 2-dimensional Brownian motion is 



conformally invariant, we can define a Brownian motion on the modular surface |2l 
Our chamber can be the standard fundamental domain; if we use the standard map 
of the modular surface to the upper half-plane, our chamber is bounded by x > —1/2, 
x < 1/2, and x 2 + y 2 > 1. Can the resulting formulas be used to compute properties 
of this Brownian motion, such as asymptotic survival probabilities, hitting times, and 
physical models? 
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